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Abstract 

We solve and characterize the Lagrange multipHers of a reaction-diffusion system 
in the Gibbs simplex of R^+^ by considering strong solutions of a system of parabolic 
variational inequalities in . Exploring properties of the two obstacles evolution prob- 
lem, we obtain and approximate a TV-system involving the characteristic functions of 
the saturated and/or degenerated phases in the nonlinear reaction terms. We also show 
continuous dependence results and we establish sufficient conditions of non-degeneracy 
for the stability of those phase subregions. 



1 Introduction 

This paper is motivated by the vector-valued reaction-diffusion equation 

dtU - AU = F{x,t,U), ing, (1) 

for U = U{x,t), defined from Q = 17 x (0, T) into M^"^^, with homogeneous Neumann 
condition on dQ x (0,T), where 17 is a bounded domain of M" and T > is arbitrary. We 
are interested in the case when every component Ui = Ui{x, t) is nonnegative and the system 
is subject to the multiphase non- voids condition with J = (1, . . . , 1) G M^"*"^: 

Af+l 

U ■ J = "^Uj = 1 in Q. (2) 

i=i 

Prom the equation ([T]) it is clear that the constraint ^ implies F{x, t,U) • J = in Q 
and so the reaction vector F should satisfy the necessary and very restrictive condition 

N N 

FN+i{x,t,V) = -^Fj{x,t,V) inQ, yV = {vi,...,Vn,l-^Vj), 0<v^<l. (3) 
i=i i=i 
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For instance, in replicator dynamics describing the evolution of certain frequencies in a 
population, one possible definition of the reaction term with this compatibility condition 
consists in choosing 

N+l 

Fi{x,t,V) =Vi[(l)i{x,t,V)-J2'"j(t>jix,t,V)] inQ, i = l,...,N + l, (4) 

j=i 

where Vi represents the i-frequency of the population and (pi the respective fitness (see, 
for instance, [TO] and [H]), the constraint ^ is essential to describe mixed strategies in 
evolutionary game theory in spatially homogeneous population dynamics (see [18] and its 
references) or to model the non- voids condition in biological tissue growing [lH [H]. In 
phase fields models, the condition ([2|) arises naturally in simulation of multiphase fiows 
(|13|) and multiphase systems with diffuse phase boundaries, as in solidification of alloys or 
in grain boundary motion (see |9] or [3]). 

Of course, in the case ([3]), in particular, if = 0, the problem becomes a simple one if 
the initial data U{0) = Uq also satisfies the constraint However the situation is entirely 
different in the general case of non trivial reactions, specially in multiphase problems where 
at least one phase "i" in a subregion of Q is absent (i.e. = 0), or fulfils another subregion 
(when Ui = 1). 

Instead of solving the system ([T]) in the Gibbs (N+l)-simplex 

Af+l 

^ = {{vi,...,VN+i) e : ^ = 1 and > 0, i = l,...,iV + 1}, 

i=i 

we shall replace this problem by the study of a unilateral problem for the vector field of the 
first N components u = (ui, . . . , un) of U, with the A'' + 1 convex constraints 

N 

Uj < 1 and Ui >0 in Q, i = I, . . . , N. (5) 

i=j 

This corresponds to solve the system of parabolic variational inequalities, at each time 

te{o,T), 

u{t) G K : / dtu{t) ■ {v - u{t)) + / Vu{t) ■ V{v - u{t)) 
Jn Jn 

> [ f{u{t))-{v-u{t)), yveK, (6) 
Jn 

under the initial condition 

tt(0) =uo = (uoi, . . . , uqn) G IK. (7) 
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Here IK denotes the convex subset of the Sobolev space H\n)^ defined by 

N 

K={ve H\n)^ : <1, Vi>0, i = l,...,N, in U}, (8) 

3=1 

where v = {vi, . . . , vn)- 

The reaction term may have a general form fi{u) = fi{x, t,U{x,t)), i = 1, . . . , with 

{x,t) G Q and U = {ui, . . . ,uj\f, 1 — ^Uj). We denote dt = — and V = i— — , . . . , — — j. 

The main part of this work is the analysis of the new unilateral problem ©-([Tl) under 
general assumptions on /: only continuity on n and integrability in (x, i) E Q. In particular, 
we prove that its solution u = u{x,t), which each component Ui satisfies a double obstacle 
problem 

0<Ui<l-^Uj in Q, i = l,...,N, (9) 

where ^^^^j denotes the sum of all — 1 components but Ui is, in fact, also the solution 
of a reaction-diffusion system in the form 

dtUi - Aui = fi{u) + /r(tx)X|^_o| 

E ^(/nH + + ing. (10) 

1< ii<. . .< ifc< JV 
i e {il, . . . , ifc} 

Here denotes the summation over all the subsets {ii, . . . , ik} of {1, ... , A^} 

1< ii<. . .< ifc< N 
! 6 {ii, . . . , ifc} 

to which i belongs, in particular, k varies from 1 to N . We also denote = g y and 
g~ = —{g AO) the positive and negative parts of a scalar function g = g~^ — g~ , Xj^ the 
characteristic function of the set A, (i.e., = I m A and Xj^ = in Q\A) and the 
characteristic function of the set 



{{x,t) G Q : {ui,+---+Ui^){x,t) = 1, Ui^{x,t) >0, j = l,...,k}, ke {!,... ,A^}. 



In particular {ui = 1} = [^{"Uj = 0}, i.e., one component is fully saturated if and only if 

the others are absent. Hence from (fTOl) we see that, in general, the respective reaction terms 
are coupled not only through the semilinear term f{u) but also through the characteristic 
functions of the saturation sets of 
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In this way, by setting for i = 1, . . . , A^, 
F,{U) = n{u) + /r(«)X|„^=o} - E ^(/nH + • • • + fM)^^H-i., 

1< ii<. . .< if.< N 
i g {il, . . . , ifc} 

N 

with U = {u,l - ^ Uj), we can solve the system ([T]) under the constraint ^ and identify 

the respective Lagrange multiphers hi = Fi{U) — fi{U) in a precise form. 

To illustrate the meaning of the system (fTQ]l . that contains 2^ — 1 + characteristic 
functions, in general, we may consider the cases = 1,2 or 3. Denoting, for simplicity, 
fi = fi{u), = we may write the Lagrange multipliers as 

h, = /rX|„^=o} - ft^l -Ufl+ h)^^12 + /3) + ^13 - Uh +h + /3) + ^123 

h2 = f2^{u,=0} - ft ^2 -Ufl+ h)^^l2 -Uf2 + /3) + ^23 " ^(/l + /2 + /3) + ^123 
h3 = fs^{u,=0} - ft^3 /3)''^13 -Uf2 + /3)+^23 " + /2 + /3) + ^123 

Ignoring the third equation and all the terms involving the third component, we may 
obtain the case N = 2. The first two terms of the right hand side of the first equation 
correspond, in the case A^ = 1, to the scalar two obstacles problem that has been proposed 
for phase separations in [H[5]. 

The mathematical treatment of this unilateral system is done in the following three 
sections. In section 2, we consider the semilinear approximation of the unique solution of 
([11)- d?!) in the case of the reaction / is in L'^{Q)^ and independent of the solution. Although 
there exists a large literature on parabolic variational inequalities (see, for instance, [16| . 
[S], [E], [I] or [8]), the direct approach of the bounded penalization used for the two 
obstacles problem in [22] (see also [IS]), extended here for the system ([TO]) , allows the use 
of monotone methods. This yields a direct way of obtaining Lewy-Stamppachia inequalities 
([26]) . obtained first by [7] for parabolic problems, implying the Wp'^ and Holder regularity 
for the solution to ((H). Similar results for the AT-membranes stationary problem have been 
obtained in [Tl|2]. We note in our case the simplification due to homogeneous Neumann 
condition. 

In section 3, we extend the existence result to general nonlinear reaction / = f{u) 
taking values in L^{Q)^ . Here we explore the fact that the convex set ^ lies in the unit 
disc and we extend the direct technique of pO]. We show also a continuous dependence 
result and, in the case of XI — f being monotone non-decreasing, in particular if / is 
Lipschitz continuous in u, also the uniqueness of solution and their strong approximation 
by the penalized solutions. 
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Finally, in the last section, we characterize the solution of the variational inequality ([6]) 
as solutions of the reaction-diffusion system (fTO]l . by extending some remarks of [23j to the 
two obstacles parabolic problem. We also show that 

k 

{ui = 0} C {Mu) < 0} and C { ^ ^ (u) > O} 

i=i 

a.e. in Q , for 1 < ii < ■ ■ ■ < ik < N , \/ k = 1, . . . , N and we can modify the system 0^ 
(see dTTj) and show that the a.e. pointwise nondegeneracy assumptions 

k 

Y^fi^{u)^0, l<ii<-- - <ik<N, k = l,...,N, 
i=i 

are sufficient conditions for the local stability of the characteristic functions ^{u-=o} 
•^iL.-ik with respect to the perturbation of the nonlinear reaction terms /. 



2 Approximation of strong solutions by semilinear problems 

In this section we consider the case where / = (/i, . . . , /jv) depends only on (x, t) and 
is given in LF'{Q)^ . 

To prove existence of solution of the variational inequality (I6])-(I71), we consider a family 
of approximating semilinear systems of equations. We define, for each e > 0, : M — > M 
by 

' if s > 

e^{s) = 1 s/e if - e < s < (11) 
-1 ifs<-e, 

and we denote 

Pu = dtu — Alt = {Pui, . . . , Pun), 

where dtU = {dtUi, . . . , dtu^) and Au = (Aui, . . . , Aujy). We also denote Pui = dtUi—Aui, 
i = 1,...,N. The approximating problems are given by the following weakly coupled 
parabolic system with Neumann condition 

Pnf + /-0,«)- Yl l{fn + + = n inQ, (12) 

1< ii<. . .< ifc< JV 
i e {ii, . . . ,ifc} 

' on5J^x(0,T), (13) 



2i 

dn 



<(0)=txoi inQ, (i = l,...,iV) (14) 
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d 

where — — is the outward normal derivative on dQ x (0, T), the meaning of > 
on ^-^ 

1< ii<. . .< N 
i S {ii, ■ ■ . ,ik} 

was explained in the introduction and 

yv = {vi,...,VN) V{ii,...,ifc} C {l,...,iV} Vij^,„i^ = Vi^ -\ \-Vi^. (15) 

Defining the penalization operator by 

N 

@,U-V=Y, [f^~Oe{u,) - (^^1 + ■ ■ ■ + - ''n...^,)]v^ (16) 

i = l 1< ii<. . .< ifc< N 

i e {»!, . . . ,ifc} 

^ 1 

j = l 1< ii<. . .< ifc< JV 

we formulate (fT2]l - (fT3]l in variational form for a.e. t £ (0,T), 

[ dtu%t)-v + [ Vu'{t)-Vv+ [ ee{u'{t))-v= [ fit)-v, yve H\n)^, (18) 
Jn Jn Jn Jn 

associated with the initial condition (fT4l) . 

Proposition 2.1. Assuming that 

f = ifi, . . . Jn) e L\Q f and uo£K, (19) 

the problem has a unique solution £ H^O,T; L^{n)'^) f] L'^{0,T; H^{n)^)). 

Proof. We begin by proving the monotonicity of the penalization operator 0£. 

In fact, recalling that 6s is monotone nondecreasing and the definition ([15]) we have 



(Belt - QeV) ■ {u-v) 

N 

1=1 



Y,f-{Oe{u^)-eM)){ 

^ ^(-^^i + •••/ij'^(6'£(l - Ui^...^k) - Ge['^-Vi^...ik)){u,^...iu -■"n.-ij, 

1< ii<. . .< z/,< A'" 

>o, 

since f~ and (/ij + • • • /i^)^ are nonnegative functions. 

The existence and uniqueness of solution E L?{^,T;H^{Q)^) is immediate by ap- 
plying the theory of monotone operators ([IE], [25])) . 
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Setting V = (uf , . . . in the approximating problem (fTSll and integrating in time, 

letting 

gl = Put = h- frOM) + E ^(^^1 + • • • + /^^''^^(l - 

1< ii<. . .< if< N 
i e {ii, . . . ,ik} 

which is bounded in L?'{Q) independently of e, we obtain that, for every < t < T, with 
= J7x (0,t), 



^ ' .e/j.\\2 I f |V7„.£|2 ^ -'^ /" l„. |2 I 1 / I £|2 , / i .e|2 



^ JO JQt ^ Jn ^ JQt ^ JQt 

The Grownwall inequality yields the uniform boundedeness (in e) of w^, first in 
L~(0,r;L2(J7)^)) and afterwards also in L'^{0,T; H\Q)^). 

Letting, formally, v = dtu^ in (fTSj) (in fact in the respective Faedo-Galerkin approxima- 
tion) and integrating in time, we get 



top 



[ \dtuf+ [ \Vu%t)\^< [ \g'\^+ [ \Vu, 
JQt Jn JQt Jn 

and so dtu^ is also bounded in L'^{Q)^ and Vm^ in L°°(0, T; L2(17)^). Therefore 

{u^}e>o isboundedin H^{Q,T;L^{Vt)^))r\L°°{Q,T;H^{Vt)^). (20) 

□ 



Proposition 2.2. Assuming (fT9|l . the solution of the problem p8|l - (|T4l) satisfies 

N 

ul>-e, i = l,...,N, ^nf<l + e. (21) 

i=l 

Proof. In fact, we are going to prove the following more general set of inequalities 

< > i = l,...,N, and uf^ < 1 + e, V 1 < ii < . . . < v < iV 

and the proof of the right hand side inequalities will be done by induction on r. 

Let us prove the case r = 1, i.e., uf < 1 + e, for all i £ {1, • ■ • Multiplying the 

i-th equation of the approximating system (fT2]l by (uf — (1 + e))~^ and integrating over 
Qt = 0, X (0, t), we have 

/ a^nf « - (1 + e))+ + / V< • V(nf - (1 + £))+ 
JQt JQt 



1< !i<. . .< N 
• 6 {ii, . . . , ifc} 
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Recalling that — 1 < < and that, in the set {uf > 1 + e}, we have 0^{uf) = and 
— n|) = —1 we get 

1 f + / |v«-(l + e))+p< / (/,_/+)(n|_(i+e))+<0, (22) 

^ JO JQ Jo 



SO (lif - (1 + e))+ = 0, i.e. uf < 1 + e. 

Assuming we have proved that u| ^ < we are going to show that u] 



We multiply the equations ij J = I, . . . ,r + I, hy - 
r + 1 and integrate over Qt. We obtain 

^ „ r+l r+l 

/ ^'<...vv+.«...vv+.-(i+^))^= / fE^.-E/^^' 



< 1+e. 



(1 + e))"^, sum from 1 to 



ll...lrlr+l 



r+l 



+E E r(/n + ---+/.j^^.(i-<...u.) «...vv..-(i+^)r 



J = l 1< ii<. . .< ifc< JV 
e {ii, . . . , ifc} 



Observe that, in the set {u[ 



ii...irir+i > 1 + e} we have uf^. > 0, for j = 1, . . . , r + 1, since, 
by induction, uf_^ = ''^fi + ' ' ' + ^t+i ~ ""fj < 1 + £• So, in that set 9£{uf.) = and, on 



the other hand, ^^(1 — 



«l...lr.«r + l ^ 



+ / |v(n; 



-1. The induction conclusion follows from 



'll...lrlr+l 



''21...?.rlr+l 



(1+^))" 



r+l 

To prove that uf > — e, we multiply the z-th equation of (fT2]l by (— — e)"*", obtaining 



- / |(-nf-e)+(t)|2+ / |V(-<-e) 



Q 



Q " 



1< ii<. . .< ifc< iV 
i e {il, . . . ,ifc} 



Let Jfc^j = {ii, . . . \ {i} and denote the elements of J^^j by j'l . . ■ jk~i- Since, in the 
set {(— n| — e)^ > 0} = {uf < —e}, we have 1 — "uf^ = 1 — n^^ — > (recall that 



<-ent)\'+ / |V(-<-£)+p < 



that implies (— u? — e)^ = 0, or uf > —e. 



-/+(-<-£)+ < 0, 
□ 
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Theorem 2.3. Assuming (fT9l) . the variational inequality JUj-jT]) has a unique solution u 
such that 

ue H\0,T;L\n)^)n L'^{0,T;H\n)^) (23) 

and 

Pu G L\Q f. (24) 

Proof. Let be the solution of the problem (fTSl) . Using the uniform estimates (in e) 
obtained in (f20]) . we know there exists u such that 



u in L (Q) strong, 



e 

vF u in L'^{0,T;H\Q)^) weak-*, 

e 

^ dtu and Pu*^ ^ Pu in L'^{Q)^ weak. 



We have it(t) G K, for a.e. t G [0,T], because satisfies the inequalities ([2T 
Given G L^(0,T;]K), set ^;(t) — w^(t) in (fTSll and integrate in time. Then 



/ 5tw^-(t;-u)+ / Vu' ■V{v-u)> [ f-iv-u'), 
Jq Jq Jq 

since / {Qeiu") - @£{v)) ■ {v - u^) < and @s{v{t)) = if v{t) G IC. Passing to the limit 

Jq 

when e — > and noting that 

lim inf / {dfu'' ■ + Vu^ ■ Vtt") > / {dfU-u + Vu- Vu) , 

Jq Jq 

we find that u satisfies ([7| and 

/ dtu-{v-u)+ / Vu-V{v-u)> / f-{v-u), G L^(0,r;lK), (25) 

Jq Jq Jq 

which is easily seen to be equivalent to jH]). The uniqueness is immediate. □ 
We remark that no regularity of the boundary dQ has been required in (fTSj) and, in 



fact, the Neumann boundary condition (I13p is only formal. In the proof of Theorem 
we have used the compactness of the sequence {w^je in L'^{Q)^ . This holds, for instance, 
for domains with Lipschitz boundaries, but also, since the sequence {u^}e is uniformly 
bounded in L°°{Q)^ , for a larger class of bounded open subsets of M^"*"^. However, the 
approximation by semilinear parabolic equations yields immediately an additional regularity 
of these strong solutions. 
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Indeed, from the definitions of 6s and G^, from (fTSjl with arbitrary f £ V{Q), > 0, 
we find 

E \{h. + --- + h,)^<P< = fi-Qe[u')<fi + f; = ft a.e. inQ. 

1< ii<. . .< if< N 
i g {il, . . . , ifc} 

(26) 

By the conclusion of Theorem 12.31 we also obtain, for each i = 1, . . . ,N, 

h- J2 \{h. + --- + fh)^<P^i<ft a.e. inQ (27) 

1< ii<. . .< ifc< JV 
i G {ii, . . . , ifc} 

and we can apply directly the second order linear parabolic theory (see [T^) in the Sobolev 

spaces 

W^'\Q) = W^^P{Q,T-LP{^)) n Lf(0,T;Ty2,P(f7))^ 1< p < oo. 
These spaces satisfy the Sobolev imbeddings, for p > (n + 2)/(2 — A;), with /c = 0, 1, 
W^^^iQ) C Ct'\Q), 0<Q<2-A:-(n + 2)/p, 

where Ca^{Q) denotes the spaces of Holder continuous functions v m. Q, with exponent a 
in the x- variables and a/2 in the i- variable and, in the case /c = 1, with Vv satisfying the 
same property (see [17], p. 80). Therefore, as a consequence of (f27l) . we conclude. □ 

Theorem 2.4. Assume that is smooth, say of class and 

f G LP{Q)^ and uq G K n VT^-^Mp^j^^JV^ 1< p < oo, (28) 

with each component uqi satisfying the compatibility condition — — = on dVL if p > 3. 

on 

Then the unique solution u of the variational inequality (l6])-([7| is such that 

u £ W^^\Qf n L°°(0,r;]K), (29) 

and, in particular, is Holder continuous in Q if p > {n + 2)/2 and has Vu also Holder 
continuous if p > n + 2. 

□ 

We observe that, when p < 2, the inclusion Wp'^{Q) C L'^{0,T; H^{Q)) only takes 
place if p > (2n + 4)/(n + 4) but, as we shall see in the next section and since K is bounded, 
(HD-dlD is solvable for any / G L^{Q)^. 
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3 Existence and uniqueness of variational solutions 

In this section, requiring the compactness of the inclusion of ff^(O) into L^(J7) by 
assuming a Lipschitz boundary 9fi, we show how we can still solve the variational inequality 
(|25l) for a more general initial condition 

N 

uoe K = {v e L^{n)^ : <1, Vi>0, i = l,...,N, in Q} (30) 

i=i 

and for general nonlinear / = f{u) defining a continuous operator from L'^{0,T;'K) in 
L^iQ)^. We shall assume that / = fix,t,v) : Q x [0, 1]^ ^ satisfies 

/ = f{x,t,v) is continuous in v for a.e. {x,t) £ Q, (31) 

3ipi£L^{Q): \f{x,t,v)\<ipi{x) V-!;G[0,1]^, for a.e. (x, t) G Q. (32) 
However, now the solution has less regularity, namely 

u G C([0,r];L2(17)^ n t) n L'^{0,T;H\n)^) (33) 

and its derivative may not be a function, since we only have 

dtu G L\Qf + L^O,T; {H\nf)'). (34) 

Hence the first term in the variational inequality (f25ll should be interpreted in the duality 
sense between L^Q)^ + L2(0, T; (ifi(O)^)') and L°°{Q)^ D L'^{0,T; H^fl)^), namely 
through the formula 

{dtu,v)t= [ Pu-v- [ Vu-Vv, yv e L'^iQ)^ n L'^{0,T;H\n)^), (35) 
jQt JQt 

for arbitrary t G (0, T] since, as we shall see, (f271) yields Pu G L}{Q)^ . 

Theorem 3.1. Under the assumptions (I30p . (I3ip and (I32p . the variational inequality (l25l) 
has a solution u satisfying ([33]) . (f34l) . (f27l) anc? ti(0) = uq we can write 

[ {Pu- f{u)) ■ {v -u) >0, V^; G L2(0,r;t). (36) 
Proof. We consider the closed convex subset of L'^{Q)^ 

N 

K = L2(0, T; K) = {-y G L^{Q)^ : > 0, i = 1, . . . , A^, J]] < 1 in Q} 

i=l 
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and we define ^ : K ^ K as the nonlinear operator that associates to each w € K the 
solution = ^{w) of the variational inequality ((25]) with / replaced hy g = f{x,t,w) 
and fixed initial data uq £ K. 

By showing that $ is a continuous and compact operator, a fixed point u = given 
by Schauder Theorem, will provide a solution with the required properties. 

Indeed, first we observe that if we consider any sequence K 3 
L'^iQ)^, by (HH) and ([32]), the Lebesgue Theorem implies 

= fiw,) > f{w) =g in L\Q)'' . 

Next, for any g G L^{Q)^ and any uq G IK we consider sequences g^ G L?'{Q)^ and 
UQy G ]K such that 



g^ > g in L^{Q)^ and uqi, > uq in L'^{^1)^ 

V 



and we denote by Uy = S{uoy,gy) the unique solution of (f25l) -(f7l) given by Theorem 12.31 
for each g^ and uq^. We observe that each component of Pu^ satisfies the inequality (f27l) 
with fi replaced by {gy)i. Prom ([251) for and u^, we easily find, for a.e. t G (0,T) 

and, integrating in time, we obtain 



sup / \u^{t) - Uy{t)\'^ + / \V{u^-Uy)\'^ < / \uQ^-uay\^ + A Ig^^-g^] 
:i<t<TJn Jo Jq Jq 



(37) 



0<t<T. 

This estimate shows that {uy}i, is a Cauchy sequence in the Banach space 

W = C([0,r];L2(O)^) n L'^{0,T;H\n)^) (38) 
with respect to the norm 

sup / \vit)\''+ [ \Vv\^y^^ (39) 
•o<t<TJn Jq ' 

and, hence, there exists a function Ug£ W 

Uy > u„ in W . 

In addition, Ug G 1?{^,T;¥.) n C7([0,r];IK) and Pug G ^^(Q)^, which implies, by 
([35]l . that 9ftig satisfies ([M]) . Hence, using ([351) . we may pass to the limit in v in 

{PUy -gy,V- Uy) = / {PUy " g y) ' {V-Uy)>0 

Jq 
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for an arbitrary v G L'^{fd,T-K) C L°°{Q)^, and using the formula 



2{dtUg,Ug)t= I \Ug{t)\^- [ \Uo\^, VtG {O , T] , 



Jn Jn 



we conclude that Ug = S{uQ,g) is the (unique) solution of the variational inequality ((251) 
(or equivalently ([36]) ) with data g S L^{Q)^ and wq £ IK. In particular, from ([37l) . we also 
obtain that, for fixed tio S K, the operator S : g i-^ iig = S{uQ,g) is Holder continuous of 
order 1/2, from ^^(Q)^ into PV". 

Since dtUg satisfies the property ([M]) . it is in fact in L\0,T;H-'{n)^), for s sufficiently 
large and, by a well known compactness embedding (see [24] or Theorem 3.11 of [25]), the 
compactness of H^{Q) C L^(J7) implies that, in fact, S regarded as an operator from 
L^{Q)^ into ii' C L'^{Q)^ is, therefore, completely continuous. Hence, $ = S o / fulfils 
the requirements of the Schauder fixed point theorem and the proof is complete. □ 

Remark 3.2. It is dear that if uq K and, in (f32l) . ipi G L?'{Q), we obtain in Theorem 
13.11 the existence of a strong solution satisfying (I23p and (l24l) . Of course, if we have the 
regularity assumptions of Theorem 12.41 i.e., G L^{Q), implying by the inequalities (l27l) 
that Pu G L'P{Q)^ , we also obtain solutions in Wp'^{Q)^ , in particular Holder continuous 
solutions if p > (n + 2)/2. 

In general (f36]l may have more than one solution, but if we assume, in addition, that 
for some A > 0, XI — f is monotone non-decreasing in [0, 1]^, i.e. 

3 A > : X\v — w - {f{x,t,v) - f{x,t,w)) ■ {v - w) > 0, ix,t) G Q, yv,w G [0,1]^, 

(40) 

in particular, if / is Lipschitz continuous in then there exists at most one solution u of 
the variational inequality (f25l) in the class ((331) and initial condition uq G IC. 

In order to prove the uniqueness of solution, we suppose that ui and U2 are two solutions 
of the variational inequality ((251) with initial condition tig G K and / = /(wi), / = f{u2) 
respectively. Then, choosing U2 and Ui as test functions, respectively, using ((401) we find 




< 



/ {f{U2) - f{Ui)) ■ {U2 -Ui) < X / \U2 - Ui 




and so, by Grownvall inequality ui = U2 a.e. in Q, since «i(0) = 1*2(0) = tto- 
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We redefine the variational formulation of the approximating problem (jlSp in the frame- 
work of this section with defined in (fTBI) and with initial condition only in L'^{Q)^ , 

[ dtu'-v+ [ Vu'-Vv+ [ ee{ul-v= [ f{u')-v, yv £ L'^{0,T;H\n)^)nL°°{Q)^. 
JQ JQ JQ JQ 

(41) 

Arguing as in Theorem l3.1l we may prove the existence of a solution of the approximating 
problem (fT2]) . with initial condition wq £ ]K as long as / satisfies (f3T]) and (f32ll . We also 
have uniqueness if we assume ([40]) . 



Theorem 3.3. Suppose that f satisfies ([3T]) . ((321) and (|40]) and uq £ K. 

Let and u be, respectively, the unique solution of the approximating problem (fT2l) 
and of the variational inequality (f25]l . both with initial condition Uq. Then there exists a 
positive constant c = c{(pi , T) such that the following estimate in the norm (l39l) of W = 
C([0,T];L2(17)^) n L\0,T;H\Q)^) holds, 

\\\u^ - u\\\ < c^/e. (42) 

Proof. We choose in (IHI) v = — u as test function. Since u G K, then 

/ ee(w^) • (H^ -U) > 

JQ 



and so 



(43) 



/ -{u^ ~u)+ [ Vu" -{u^ -u)< [ f{u') ■ {u" - u). 

JQt JQt JQt 

Choosing, as test function in ([25]) = ((uf — • • • , (u^ — we get 

/ dtu • (li^ - u) + / Vit • V(u^ - It) 
JQt JQt 

> I f{u) ■{u'-u)+ [ [Pu - f{u)] ■ {u' - v') (44) 
JQt JQt 



and subtracting (|44|1 from (|43]) we get 

^ Jn JQt 



<[ {f{u')-f{u))-iu'-u)+f [Pu-f{u)]-{v'-u') (45) 
JQt JQt 



<X |u^-ixp + e / \Pu-f{u) 
JQt JQt 
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since — u^\\iaa(^Q-^N < e. Letting C = C{ipi,T) = \\Pu — f {u)\\ l'^ [q)n and dropping 

the nonnegative term / |V(w'^ ~ '^)\'^ in ([45]) we obtain, by application of the Grownwall 

jQt 

inequahty, 

„2At 



/ \u'{t) -u{t)\'^ < 2eCe^ 
Jn 



and using again (I45]l . also 

\\\u^ — u\\\ < cxfe. 



□ 



With similar arguments we may give a continuous dependence result for solutions of the 
variational inequality (l36ll . 

Suppose we have a sequence > / in the following sense 

V 

_ jv(x^t^v) are continuous in i> G [0,1]^, for a.e. G Q 1 

f{:r,v) — /(•,-,^;) in Li(Q)^ for all fixed 17 G [0,1]^. J ^^^^ 

In addition, the assumption (l32]l is satisfied for all / uniformly in i.e., there is a 
common such that ([32]) holds for all and the initial data are such that 

^3ul > -uo in L^{^)^ . (47) 

Hence, by Theorem 13. 11 it is clear that there are solutions {w'^j^gN to the corresponding 
problems associated with and and, moreover, they satisfy ([33]) and ([34]) uniformly in 
z/, i.e., their norms in those spaces are bounded by a constant independent of v. Therefore, 
we have a function u in the same class ([33]) and ([M]) . and a subsequence, still denoted by 
z/, such that 



vy u in L2(0,r;i/i(O)^) weak and in L°°(0,r;K) weak-* (48) 

u" >u a.e. in Q and in L'P{Q)^ , Vl<p<oo. (49) 

V 

By assumption (j46]) and Lebesgue Theorem, we conclude first that f^iu") > fi'^^ 

V 

a.e. in Q and in L^{Q)^ , as well as 

f f-^u^ ■ u'' > [ fiu) ■ u, (50) 

Jq ^ Jq 

[ {r{un - f{u)) . {u^ - u) — > 0, (51) 

JQ 
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since, in particular, \u'^\ < 1 and \u\ < 1 a.e. in Q. 
Recalling (f27l) for each i/, we may take the limit in 

/ {Pu'' - r{u^)) ■{v-u'')>Q (52) 
JQ 

for a fixed v G L^(0, T;]fC). Using ([50]) and (l48ll . that in particular imply 

PueL^iQ)'^ and liminf / Pu" ■ u" > I Pu ■ u, VtG(0,r), 



we conclude that u is a solution of (|36|1 with initial condition Uq- 

Using V = w^(o,t) + u^^(t,T) in ([52]l and t; = u'^^(o.t) + '^^(t,T) ™ pSil we find, for a.e. 

t G (o,r), 

and, by ([5T]) . we conclude that u'^ > u strongly in W . Therefore, we have proved the 

following result 



Theorem 3.4. Ifu^ denotes the solution to the variational inequality (|36p with satisfying 
the assumptions (I46p and (I32p uniformly in v and initial condition satisfying (|47p . i/ien ^/lere 
ea;isfe a subsequence {u'^}ueN such that 

>u m C{[0,T];L^{n)^ nK)nL'^{0,T;H^{n)^)nLP{Q)^, Vl<p<oo, 

where u is a solution to (|36p corresponding to the limit f and the limit initial condition Uq. 
In addition, if f satisfies (l40l) . by uniqueness of u, the whole sequence {u'^}^(zfq converges. 

□ 



4 The multiphases system and its characterization 

In this section we consider a variational solution u of ([251) obtained in Theorem l3.ll i.e., 
satisfying ((331) and (fMp . Setting 

Wi{u) = l-^Ui, i = l,...,N, (53) 
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each component Ui satisfies a double obstacle problem 

< Ui{x,t) < Wi{x,t) a.e.{x,t)eQ, i = l,...,N. (54) 

For an arbitrary nonnegative and bounded function if = (p{x,t) defined for {x,t) G Q, 
such that 

K'^ = {ve L'^{0,T;H^{Q) : < v < if in Q} / 0, (55) 
and for a given g £ L^{Q), we may introduce the parabolic double obstacle scalar problem 

u e IfC^ : / dtu{v -u)+ I Vu-V{v-u)> / g{v - u) e K^, (56) 

Jq Jq Jq 

subject to a given compatible initial condition 

u{0) = no in 0. (57) 

For each i = I, . . . , N , we have Ui G ICq ' and, by choosing in (f25l) v G L^(0, T; IK), such 
that Vj = Uj for j ^ i and Vi = v £ Kq' arbitrarily, it is clear that Ui is a solution of 
the scalar double obstacle problem ([56]) with ip = Wi and g = fi{u). Hence we can obtain 
further properties of our solution by applying the general theory of the obstacle problem. 
For the sake of completeness we prove here the result below. 

Let 

^eL^{id,T;H^{n))r\L°°{Q) with > a.e. in Q, (58) 
dt^&L^{{),T]{H^{n))') with Pip e L\Q), |^ = on x (0, T), (59) 

and 

g G L^iQ), uo G L^{^), < uq < fiO) in n. (60) 
We observe that (f59l) means that if satisfies the formula 

{dtip,v)t= [ vP^- [ V(p-Vv, yv G L'^{o,T-H\n))nL°^{Q). 

jQt JQt 

Proposition 4.1. Under the assumptions (j58p - (|60|) the unique solution u G Kq to the scalar 
problem ([56]) - ([57l) is such that 

u£ C{[0,T];L\n)) n L°°(Q), dtu £ L\Q) + L\0,T- {H\ny) , (61) 

and it satisfies the parabolic semilinear equation 

Pu = g + g'^{u=Q} - {P^ - g)'^{u=^} a.e.inQ. (62) 
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Proof. Using the function 6s given by (fTTI) and defining 

Mv) = 9-0e{v) - {P^ - g)~ee{y^ - v) (63) 
we can consider the approximating problem, for e > 0, 

/ {Pu" + ^eiu'))v= [ gv, VvG L^{0,T;H\n)) n L°°{Q), (64) 
Jq Jq 

with the initial condition n^(0) = uq in $7. Since •&£ is monotone and ip is bounded, arguing 
as in Theorem 13.11 the problem ([Ml) has a unique solution u'^ in the class (f6T]) . Moreover, 
it satisfies 

— e<u'^<ip + e a.e. in Q, (65) 

as we can show by choosing, in (f64ll . v = {—u^ — e)~^ and v = {u^ — ip — e)+, respectively. 
Indeed, in the first case we have 



vPv = - / vPu"" = / V (??(u^) -g)= {-g~ -g)<0, 

Q Jq J{v>0} J{u'^<-e} 

since ■i?e(n'^) = — 1 and i^e^^p — Ug) = 0, because < —e and (p — u^ > e, and, in the second 
case. 



/ vPv = / vPiu"" -(p) = [ v{g- i?(ii^) - Pip) 
Jq Jq J{v>f)} 

= [ {-{P^ - g) - {P^ - g)-) < 0, 



' {ip—u^>e} 

since "dei^ — W^) = —1 and 'i?e(u^) = if ip — Ue < —e and > (p + e. 

Hence, using the monotonicity argument, we easily conclude that u = liniu'^ G Kq is 
the unique solution of the variational inequality (f56]l . Remarking that, from ([631) we have 

-g~ < i?e(^i^) < {P'P - g)~ a.e. in Q, 
from ((Ml) we deduce in the limit the Lewy-Stampacchia inequalities 

(P(/j — g)^ < Pu — g < g^ a.e. in Q. 
In particular, this yields Pu G L^{Q) and ([56]) implies that u also solves 

[Pu- g)[v -u)>Q, VfGt^, (66) 

where = {w E L2(Q) : < v < 99 in Q} C L°°((5). 
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Let O C Q be an arbitrary measurable set and set v = u in Q \ and v = 6ip in O , 
with 6 £ [0, 1], in (f66l) . Since O is arbitrary, we conclude the pointwise inequality 

{Pu-g){(l)-u)>0 y(t)e[0,ip{x,t)] a.e. in Q, (67) 

which implies, up to null measure subsets of Q, 

Pu-g>0 in {u = 0}, Pu - g < in {u = ip}, (68) 

Pu = g in A = {0 < u < ip}. (69) 

On the other hand, arguing as in Lemma 2 of [23] and noting that V = (u, —Vu) G 
Li(Q)"+i and D-V = Pue L\Q), with D = {dt,d^,, . . . , d^J, we have 

Pu = a.e. in {u = 0} and Pu = Pep a.e. in {u = if}. 

Hence, by ([67l) . up to neglectable sets, we have {u = 0} C {g < 0} and {u = ip} C 
{Pf < g}, and using also ([681) . we finally conclude ((621) . □ 

Theorem 4.2. Any solutions u of the variational inequality (I25p (or (l36l) ) under the con- 
ditions of Theorem 13.11 satisfy the semilinear parabolic system 

Pui = fi{u) + /r(M)X{„_o| 

E \U^An) + ■■■ + f^,{'^)r^n...^, a.e. tn (70) 

1< ii<. . .< ij.< JV 
» 6 {ii, . . . , ifc} 

where = ^/^^ , for k = 1, . . . , N , denotes the characteristic function of 

Iii...ik = {{x,t) G Q : Ui^,„i^{x,t) = 1, Ui^{x,t) > for all j = I, . . . , k}. (71) 

Proof. We notice that the regularity (f58ll . ([591) . holds for if, = 1 — ^^tij, so can be 

chosen as the upper obstacle of each component Ui, i = 1, . . . , N , of u, to which we can 
apply the conclusions of Proposition 14.11 Since {ui = 0} C {/i(w) < 0} a.e., for each 
i = 1, . . . , A^, we have 

Pui = fiiu) + fi{u)X^^_Q^ - {Pwi - /j(w))"^{«,=^,,„,>o} in Q, (72) 

and the condition ([70]) will follow if we show that 

(Pwi - /i('a))"X|„_^. = i(/.^(t^) + ... + /.^(^))+Xi^ in Q, 

1< ii<. . .< ifS N 
i g {ii, . . . ,ifc} 

(73) 
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Observe that 

{ui = Wi, Ui > 0} = IJ Ih...ik, 

1< zi<. . .< if< N 
i e {ii, . . . ,ifc} 

and these sets are a.e. disjoint. Here the union is taken also over all the subsets {ii, . . . , ik} 
of {1, ... , A^} that include i and over all A; = 1, . . . , A^. We remark that Pwi = Pui in that 
subset and 

• in the sets Ii = {ui = 1}, Pwi = and [Pwi — fi(u))~ = fi{u)^, for i = 1, . . . , N; 

• in each set for > 2, as we shall see, 

(Pui - h{u))- = \{fM + ■■■ + fi,{u)t^ 

and this fact concludes the proof. 

Let {xQ^to) G Iii...ik- Recall that {ii, . . . ,ik} is the set of indexes for which we have 
< Ui^ (xo, to) (notice that i £ {ii, . . . , ik})- Denoting a = min{nj^. (xq, to) : J = 1; • • • > k}, 
the set O = f]^^i{ui^ > a/2} is measurable and contains {xo,to). Given any measurable 
set u; C O, choose, in ([36]) . as test function v = {vi, . . . ,vn) defined by 

■Uii = tiii ± (5X^, = Ui. =F SX^ for a fixed j e {2,. . . ,k}, vi = ui'il ^ 

observing that 

N N N 

Vj = Uj ± 5X^ =p SX^ = Uj < 1 
i=i i=i j=i 

and 

Vj >(), j = I, . . . , N, as long as < < a/2. 
Returning to the inequality (f36l) and setting Sj = Puj — fj{u)^ we get 

±S [ Si,X^T6 [ S^^X^>0. 
JQ JQ 

Since to D {(xc^o)} was taken arbitrarily in O and (xo,to) is a generic point of Ii^...ik^ 
we conclude that 

Si^=Si., s,.e. m. Ii^,„ik, for any j e {2, . . . , fc}. (74) 

N 

Recalling that ^^Puj = Pui-^„,i^, = 0, in the set Ii-^...i^ we get, using ([741) . that 
kSi^ = Si^^ VSi^ = {Pui^-fi^)-\ \-{Pui^-fi^) = Pui-\ \-PuN-{fii-\ h/ifc), 
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where, for simplicity, we set fj = fj{u), and so 

Si = Si^ = ~'^{fii + ■ ■ ■ + /ifc) • 

But in we have 5^ < (recall that Ui = Wi and (f68l) ) and so 

(Pui - fi{u))- = -{Pui - fi{u)) = -5i = ^ (/,, + ••• + /,J = + . . . + + . 

□ 

Corollary 4.3. Let u be the solution of the variational inequality ((25]) (or ((36]) ) under the 
conditions of Theorem 13. 1[ 

Then, denoting by \A\ the (n + 1)-Lebesgue measure of Ac Q, we have 
k k 
\{Y.f^,{u) <0} n {Y,U^, = l,U,^>0,j = l,...,k}\=0 (75) 

/or eac/i partial coincidence subset hx.-.i,,, as well as 

|{/,(u) >0}n{^z, = 0}| =0, i = l,...,N. (76) 

Proof. Being defined in ((TT]) . using the equation ((TQl) . we obtain, for each with 

j = l,...,k, denoting fi^ = fi^ (u), 

Puij = fi^ -^{fii^ ^ fik)~^ a.e in 

Summing these k equations, we have 

k 

= = ^ - (/ii H K /ifc)"^ = (/ii H ^ fik)~ a.e in 

A; 

So, in = { Z]j=i = 1) ""jj > 0, j = 1, . . . , A;} we have ^ /j^ > a.e. and ({TSj) 

i=i 

follows. 

The proof of ((76]) is similar (recall ((68l) ). 

□ 



As a consequence of this corollary the semilinear system ((70() can, in fact, be written in 
the equivalent form for i = 1, . . . , N , 

Pui = fi{u) - fi{u)X^^_Qy 

E lifni'^) + ■■■ + f^d^))^h■■■^. a.e. inQ, (77) 

1< ii<. . .< if.< N 
i g {il, . . . , ifc} 
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k 

since {ui = 0} C {fi{u) < 0} and Iii...ik { fiji'^) ^ 0} up to a neglectable subset of 

Q. 

This remark combined with the continuous dependence of the variational solutions ob- 
tained in Theorem 13.41 yields an interesting criteria of local stability of the characteristic 
functions of the coincidence sets in the Lebesgue measure. Denote 

= ''^{<^...ij^=l, ■"r>0 Vj=l,...,fe}) I <ii < ■ ■ ■ < ik < N, i e {ii, . . . 

Theorem 4.4. Let the assumptions and notations of Theorem 13.41 hold. Suppose that in 
some subset of positive measure lo Q Q the following assumption on the limit problem holds 

k 

^/^^(u)/0 a.e. inuj, I < ii < ■ ■ ■ < ik < N, k = l,...,N. (78) 
i=i 

Then the associated characteristic functions are such that 

X|„.^o} ^K=o} inLP{Lo), Vi = l,...,iV, (79) 

inU'iio), Vii,...,4, (80) 

for all p, 1 < p < oo. 

Proof. We observe that each u'^ solves the system 

P< = /r - /r^K=o} - E l(fn+--- + fO K..^. a.e. in Q (81) 

1< ii<. . .< ik< N 
i G {ii, . . . , ifc} 

where, for simplicity, we set f^ = f^{u^). By the convergence > u, we have 

Pu"^ ^ Pu in the distributional sense. Since < ^ij^,„i^, ^ 1) there exists with 

< < 1 in Q, such that 

in L°-iQ) weak- * . 
Analogously, for some with < X*q < 1 in Q, we have 

^{«^=o} — ^ K,o in L°°{Q) weak- * . 
We are going to prove that, in fact, 

Kfl = ^{u,=o} and X*^ = a.e. in w. 
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which concludes the proof, since the weak convergence to characteristic functions in I^{u)) 
is in fact strong, as it is well known. 

Passing to the limit in (fST]) . we obtain 

Pui = fi - fiKfl - ^{fii + --- + fik)^ii...h a.e. in Q 

1< ii<. . .< ifc< JV 
i e {ii, . . . 

where, for simplicity, we have also set fi^ = fi^{u). 

But each Ui also solves the equation ([77l) . so, by subtraction, we obtain a.e. in Q, 

- /.(^K=o} - ^:,o) - E l(fn +■■■ + - = 0- (82) 

1< ii<. . .< JV 
i 6 {il, . . . , ifc} 

Noticing that Xj^i/^QjuJ' = 0, passing to the limit, we get ^*fiUi = 0, which means that 
X*o = whenever Ui > 0. To conclude that X*q = Xju.=o} we only need to prove that 
X*Q = 1 if tij = 0. 

Recall that the sets {uj = 0} and Iii...ik^ 1 < < • • • < ^fc < A^, i G {ii, . . . ,ik}, 
k = 1, . . . , N , are mutually disjoint. Hence in {ui = 0} we obtain 

-/.(I - ^:,o) + E ^-^^1 + ■ ■ ■ + = 

1< ii<. . .< z/,< A'" 
i 6 {ii, . . . , ifc} 

and since the left hand side is nonnegative, by the assumption ([781) we conclude that 

X*o = l and K...i,=0 inK = 0}nu;. 

Since — ~ ^ a.e. in Q, taking the limit in z^, we also obtain 

'^ii...ife(l ~ = a.e in Q, i.e. = if Ui-^,,,i^. < 1. It remains to evaluate 

'^ii...ik when Ui^,,,ik = 1 and Ui^ > 0, for all j = 1,...,A; or when Ui. = 0, for some 
j = l,...,k. 

In this later case, where Ui^ = 0, for some j = 1, . . . , A;, we have ^ii...ik = and, since 
we already know that Xj^^ =o} = ^ij ,o j from ([82]l for the index ij , we get 

1< ii<. . .< ifc< iV 
g {il, . . . , ifc} 

Then, by the assumption ([78]) we have = in [Q \ Iii...ik) ^ ^■ 
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Finally, in again from ([82]) . we obtain 

^(/n + --- + /.J(l-^n...J=0 
and the assumption ([78]l yields that = 1, completing the proof. □ 
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